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ABSTRACT 

We present near-infrared emission-line flux distributions, excitation and kinematics, as 
well as stellar kinematics, of the inner 520 x 520 pc 2 of the Seyfert 2 galaxy NGC 5929. The 
observations were performed with the Gemini’s Near-Infrared Integral Field Spectrograph 
(NIFS) at a spatial resolution of - 20 pc and spectral resolution of 40 kms' 1 in the J- and 
Ki-bands. The flux distributions of H 2 , [Fen], [Pn], and H recombination lines are extended 
over most of the field of view, with the highest intensity levels observed along PA = 60/240°, 
and well correlated with the radio emission. The H 2 and [Fe 11 ] line emission are originated in 
thermal processes, mainly due to heating of the gas by X-rays from the central Active Galactic 
Nucleus (AGN). Contribution of shocks due to the radio jet is observed at locations co-spatial 
with the radio hotspots at 0750 northeast and 0760 southwest of the nucleus, as evidenced 
by the emission-line ratio and gas kinematics. The stellar kinematics shows rotation with an 
amplitude at 250pc from the nucleus of -200 km s -1 after corrected for the inferred inclination 
of 18.3°. The stellar velocity dispersion obtained from the integrated K-band spectrum is <x* = 
133 + 8 km s -1 , which implying on a mass for the supermassive black hole of M. = 5-2-L2 x 
10 7 M 0 , using the M. - <x* relation. The gas kinematics present three components: (1) gas in 
the plane of the galaxy in counter-rotation relative to the stars; (2) an outflow perpendicular 
to the radio jet that seems to be due to an equatorial AGN outflow; (3) turbulence of the 
gas observed in association with the radio hot spots, supporting an interaction of the radio 
jet with the gas of the disk. We estimated the mass of ionized and warm molecular gas of 
- 1.3 x 10 6 M 0 and -470 M 0 , respectively. 

Key words: galaxies: individual (NGC 5929) - galaxies: Seyfert - galaxies: ISM - infrared: 
galaxies 


1 INTRODUCTION 

Detailed mapping of the gas distribution, excitation and kine¬ 
matics around Active Galactic Nuclei (AGN) is fundamental 
to the understanding of the physics behind the AGN feeding 
and feedback processes. Near-IR integral field spectroscopy 
(IFS) on 8-10 m telescopes of nearby galaxies has become a 
powerful tool to better understand these processes, as it pro¬ 
vides a two-dimensional coverage with spatial resolution of 
a few to tens of parsecs at a spectral resolution that allows 
to resolve the gas kinematics and permitting isolating inflows 
and outflows in the central kiloparsec, where the relevant 


processes occur fe 
2009; Davies et a 


Kiireieiai. zuus; iviuiier-^ancnez et at. 

2009; Storchi-Bergmann et al. 201C; 

Schartmann et al. 

20101: Riffel & Storchi-Bergmann 2011a; 


Riffe l. Storchi-Bergmann & Winge 20131 : iDavies et~afl l20l4 
Riffel. Storch i-Bergmann & Riffel 20 14bl) . 


Some of the main results obtained from IFS of the inner kpc 
of active galaxies are: (I) ionized and molecular gas have distinct 
flux distributions and kinematics. The near-IR line-emission at 
these scales is originated by the heating and ionization of the 
gas by the AGN radiation and shocks produced by radio jets 
(e.g. IRiffel et al. 20068: iRiffel. Storchi-Bergmann & Nasarl 1201 Oal : 
iRiffel. Storchi-Bergmann & Wingel 2013h . The molecular H 2 
emitting gas is usually more restricted to the plane of galaxies, 
with kinematics dominated by rotation and inflows in the disk 
in most cases. The ionized gas emission traces a more disturbed 
medium, usually associated to outflows from the AGN, but 
frequently showing also a rotation component from the disk 
of the galaxy (e.g. iRiff el. St orc hi-Ber gmann & Naga r 2010a ; 


Riffel & S torchi-Be rgmannl 2011 alibi: IStorchi-Bergmann et al.l 

2010; Iserlohe et al. 2013; Mazzalay et al. 2014). (II) Outflows 
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are seen within hollow cones or from compact structures with 
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mass outflow rates in the range 10 -2 - 10 1 M 0 yr _1 for low- 
luminosity Seyfert and Low-Ioni zation Nuclear Emission-line 
Regio ns (LINERs) galaxies (e.g. iRiffel et al.ll20()3 : ISchonel et al.1 
2014) and 10 2 - 10 3 M©yr _1 for high-luminosity Seyferts (e.g. 


McElrov et al. 1 120141) . The outflows are obser ved with velocities 


from 200 km s 


_ _ to 800 km s l in ionized gas (iStorchi-Bergmannl 

l2014ah and -150 km s" 1 in H 2 when present ( Davies et al. 2 0141) . 
(Ill) Inflows are observed in H 2 (e.g. iMuller- Sanchez et al. 2006; 


IRiffel. Storchi-Bergmann & Wingell2013l : lMazzalav et al.l2014l) in 

Seyfe r t galaxies and in low-ion i zation gas in LINERs (IFathi et al.l 

20061: IStorchi-Bergmann et akl 12007 : Schnorr-Miiller al] l201ll. 


2014 ^13) with mass inflow rates in the range 10 -1 - 10 1 M©yr _1 
Jstorchi-Bergmannl201 4al lbh . (IV) The stellar kinematics in Seyfert 
galaxies reveal cold nuclear structures composed of young stars, 
usual l y associated with a significant gas reservoir (iHicks et al.l 


120131: [Riffel et al1l2010cl. 1201 id: IStorchi-Bergmann et al.ll2012f)~ 

NGC 5929 is a spiral galaxy with a Seyfert 2 nucleus located 
at a distance of 35.9 Mpcfl. NGC 5929 has a companion, NGC 5930 
at 20" northeast from it ( Pagdl 19521) . It presents a well defined bi¬ 
polar radio jet oriented along the position angle PA^60°, showing a 
triple structure with two bright hot spots, one located at 075 north¬ 
east from the nucleus and the other at 076 south-west from it. The 
third an d fainter radio structure is observed at the nucleus of the 
galax y dUlvestad & WilsoiJ Il984l: Iwilson & KeelIl989l: ISu et al.1 
ll996h . The orientation of the m ajor axis of the large scale disk is 
PA^45°(e.g. lSchmitt et al.lll997l) . 

The gas kinematics of NGC 5929 has been studied for at least 
3 decades. Using long slit spectroscopy along the orientation of 
the radio jet, iKeell dl985h concluded that the ionized gas kine¬ 
matics is consisten with rotation in the galaxy disk. [Whittle et al.1 
( 1986 ) obtained long-slit spectra along PA=60° and PA=-30° cov¬ 
ering the RJ3 and 0 111 ] emission lines and found double-component 
emission-line profiles along PA=-30°. They suggested that the ori¬ 
gin of these profiles was the superposition of the two components 
observed along the PA=60° (one in blueshift to the north-east and 
one in redshift to the sout h-west) and unresolved by their observa¬ 
tions, at a seeing of 176. [Rosario et al.1 (1201 Ol) report the detection of 
shocked gas associated to the north-east radio hot spot by compar¬ 
ing optical long slit spectra obtained with the Hubble Space Tele¬ 
scope Imaging Spectrograph (HST/STIS) with radio images. These 
authors found that the low ionization gas (traced by HJ3 emission) 
shows broader emission-line profiles than the high ionization gas 
(traced by [O iii]/15007 emission line) at locations near the radio 
hot spot and suggested that this broadening is due to the interac¬ 
tion of the radio jet with the NLR. Evidence of shocks is also ob¬ 
served at the location of the south-western radio hot spot, for which 
[Ferruit etld . (1999 ) showed that shock models with velocities of 
~300 km s 1 are able to reproduce the optical and UV emission-line 
ratios. The low spectral and spatial resolutions and/or small spatial 
coverage of the data have not allowed the studies above to reveal 
the complete scenario for the gas kinematics in the nuclear region 
of NGC 5929, what we have been able to do now with our NIFS 
observati ons. _ 

In Riffel, Storchi-Bergmann & Riffejj (|20 1 4bl , hereafter called 
Paper I) we have reported already the discovery of a peculiar struc¬ 
ture along PA=-30/150° in this galaxy from these observations, 
where the emission lines present double components. This structure 
was interpreted as being due to the interaction of the acquired gas 


1 as quoted in NASA/IPAC EXTRAGALACTIC DATABASE - 
http ://ned. ipac. caltech. edu/ 


from the companion with an “equatorial outflow” from the accre¬ 
tion disk around the central supermassive black hole (SMBH). In 
the present paper, we present a more complete study of the overall 
gas kinematics, as well as of the gas excitation and stellar kinemat¬ 
ics of the inner 520 x 520 pc 2 of NGC 5929. 

NGC 5929 represents the tenth galaxy of a series ob¬ 
served by our group AGNIFS (AGN Integral Field Spec¬ 
troscopy) in which we present the emission-line flux distri¬ 
butions and kinematics as well as the stellar kinematics of 
the inner few hundreds of parsec of nearby active galaxies. 
The previous papers of this series have presented similar jrtud- 
iesof the following ga laxies: I - ESQ4 28-G14 (IRiffel et al.l 


(Rdfel et al. 

2009): IV - NGC 4151 (fstorchi-Bergmann et al. 

2 O 0 H 2 OIO 

); V - Mrk 1066 IRiffel, Storchi-Bergmann & Nagar 

2010a|;IRiffel & Storchi-Bergmann 201 la; Riffel et al]20103); VI - 

Mrk 1157 (Riffel & Storchi-Bergmann 2011J Riffel et al. 20113); 


YII - Mrk79 (IRiffel. Storchi-Bergmann & Wingd 120131) : VIII 
- NGC 1068 (IRiffel et al.l l2014al: IStorchi-Bergmann etail 120121: 
iBarbosa et al.ll20l4) : IX - Mrk 766 ( Schonel et al .120 1 T 

This paper is organized as follows: In Section 2 we present 
the description of the observations and data reduction. In Sec. 3 we 
present the results for the continuum and line emission, including 
line-ratio maps as well as the gas and stellar kinematics. The dis¬ 
cussion of the results is presented in Sec. 4 and the conclusions are 
presented in Sec. 5. 


2 OBSERVATIONS AND DATA REDUCTION 

NGC 5929 was observed using _the Gemini Near-infrared Integral 
Field Spectrograph (NIFS - iMcGregor et al.ll2003l) on the Gemini 
North telescope operating with the adaptive optics module ALTAIR 
on the nights of March 16, May 24 and June 16, 2011 under the ob¬ 
servation Programme GN-2011A-Q-43. The observations covered 
the J and Ki spectral bands, resulting in a wavelength coverage from 
1.14yum to 136pm and 2.10pm to 2.54^m, respectively. The total 
on source exposure time for each band was 6000 s and the obser¬ 
vations were splited into 10 individual on source exposures plus 5 
sky exposures for each band. 

The data reduction followed standard procedures and was ac¬ 
complished using tasks contained in the nifs.gemini package which 
is part of iraf software. The procedures included the trimming 
of the images, flat-fielding, sky subtraction, wavelength and s- 
distortion calibrations, remotion of the telluric absorptions and flux 
calibration by interpolating a black body function to the spectrum 
of the telluric standard star. 

The final data cubes contain ~4500 spectra at a angular sam¬ 
pling of 0705 x0705, covering the inner 3"x3" ( 520x520 pc 2 ) of 
NGC 5929. The angular resolution is 0712 for both bands (corre¬ 
sponding to ~20 pc at the galaxy), as obtained from the full width at 
half maximum (FWHM) of the images of the telluric standard stars. 
From the FWHM of the emission-line profiles of the Ar lamps we 
obtain the spectral resolution of 1.7 A for the J band and 3.2 A for 
the K band, corresponding to velocity resolutions of ~40 km s -1 for 
both bands. 

In the top-left panel of Fig.Q] we present a large scale R- 
band image of NGC 5929 and its companion NGC 5930 from the 
Sloan Digital Sky Survey, covering 1.5'xl. 5'(Baiflard e t al. 1120111) . 
The top-right panel shows an optical image of NGC 5929 ob¬ 
tained with the Hubble Space Telescope (HST) Wide Field Plan¬ 
etary Camera 2 (WFPC2) through the broad-band filter F606W 
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Figure 1 . Top left: Large scale image of NGC 5929 and NGC 5930 obtained in the R-band from the SDSS DR4jBaillard e tap201ll) . Top right: Optical image 
of the inner 30"x30"of NGC 5929 obtained with the HST WFPC2 trhought the filter F606W iMalkan. Goriian <^Rnjn^9ij) . with the NIFS field of view 
represented by the green box. Middle: NIFS J-band sample spectra for the central region of NGC 5929, with the main emission lines identified. Bottom: NIFS 
K-band sample spectra for the central region of NGC 5929 with the main emission/absorption lines identified. The spectra shown were integrated within an 
aperture of 0725x0725 and the positions selected correspond to the south-western radio hotspot (0760 SW of the nucleus), for the nucleus (NUC), the position 
of the north-eastern radio hotspot (0750 NE) and for two positions along PA=-30° (0735 NW and 0735 SE), along of the peculiar outflow reported in Paper I. 
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from lMalkan. Goriian & Tam (1998), covering the inner 30"x30". 
The green box represents the field of view of our NIFS data. 
The middle panel shows sample spectra obtained from the NIFS 
J-band datacube within apertures of 0705x0705 for the follow¬ 
ing positions (from top to bottom): 076 southwest of the nu¬ 
cleus, corresp onding to the position of the southwestern ra¬ 
dio hotspot (lUlvestad & Wilsonl 1 19891) : 075 northeast of the nu¬ 
cleus, c orresponding to the posi tion of the northeastern radio 
hotspot (lUlvestad & Wilsonl 1 19891) : the nucleus; 0735 northwest 
and 0735 southeast of the nucleus along PA = -30°, perpen¬ 
dicular to the radio jet, where we found that the emiss ion lines 
presen t double components telffeL Storchi-Bergmann & Riff ell 
l2014bl) . The following emission lines are identified in the J-band 
spectra: [P n]41.14713, [P n]41. 18861 yum, [Fen]41.25702/im, 
[Fen]41.27069 yum, [Feii]41.27912yum, Hi Pa£T1.28216/im, 
[Fen] 41.29462 yi/m and [Fen]41.32092 fim. The bottom panel of 
Fig. in shows the spectra IN the K-band at the same positions of 
those IN the J-band. We identified the Hi Bry42.16612 jim emis¬ 
sion line and the H 2 emission lines at 2.12183, 2.15420, 2.22344, 
2.24776, 2.40847, 2.41367, 2.42180, 2.43697 and 2.45485 //m. The 
CO absorption band-heads at ~2.3 yum used to measure the stellar 
kinematics are also identified in the K-band spectra. 
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Figure 2. Sample of the [Fen]/U.2570, H 2 42.1218 and Pa/? emission-line 
profiles at three locations: 0735 southeast of the nucleus, along the SE-NW 
strip at PA = -30° (top panel); 0750 northeast of the nucleus, at the position 
of a radio hot spot (middle panel) and 0760 southwest of the nucleus, at the 
position of the other radio hot spot (bottom panel). The Pa/? line profile 
is shown as a continuous line, the H 2 as a dotted line and the dashed line 
represents the [Fe ii] emission line profile. 


3 RESULTS 

3.1 Emission-line profiles 

As discussed in Paper I, the emission lines present double com¬ 
ponents along a strip (hereafter identified as the SE-NW strip) 
crossing the nucleus at PA = -30°, which have been attributed 
to gas outflowing from the nucleus perpendicularly to the radio 
jet. In Figure |2] we show a sample of the emission-line profiles of 
[Fen]/11.2570, H 2 42.1218 and Pa/? from three locations: along the 
SE-NW strip at 0735 southeast of the nucleus (top panel) and at 
the positions of the radio knots along PA = 60° at 0750 north-east 
(middle panel) and at 0760 south-west (bottom panel) of the nu¬ 
cleus. The double components, reported in Paper I, are evident at 
the position 0735 south-east, not only for [Fen] but also for the 
other emission lines, with profiles presenting similar widths. In¬ 
deed, all emission lines detected along the SE-NW strip present 
similar double components. At the positions of the radio hotspots or 
knots, the [Fe n] profile is the broadest, H 2 is the narrowest and the 
Pa/?one is in between the two.The H 2 and Pa/? are well reproduced 
by a single Gaussian curve, while the [Fe n] profile clearly shows 
the presence of more than one component. A blue wing in the [Fe n] 
profile is observed for the position at 07 50 northeast of the nucleus 
and at 0760 southwest of the nucleus its profile is broader than a 
Gaussian curve. At positions away from the strip at PA = -30° and 
distant from the radio knots the emission lines are well reproduced 
by a single Gaussian component, as discussed in Paper I. 

As shown in Fig. \T\ several emission lines are observed in 
the J- and K-band spectra of NGC 5929. In order to measure the 
emission line flux distributions we fitted the observed line profiles 
at e ach sp atial position by Gaussian curves using the profit rou¬ 
tine |RjjM ( 2010b ). In table[l]we present the measured fluxes for 
the emission lines within an aperture of 0735 x0735 at five posi¬ 
tions: the locations of the radio hotspots, at 0750 northeast and 
0760 southwest of the nucleus; the nucleus; at 0735 southeast and 
0735 northwest of the nucleus along the strip at PA = -30° (re¬ 
gion where the equatorial outflow was observed in Paper I). Val¬ 
ues followed by were obtained by the sum of the fluxes of the 


two components. The [Fen]41.2570 line profile was fitted by three 
Gaussian curves at 0760 southeast of the nucleus (at the location of 
a radio hotspot), and the flux quoted in the table is the sum of the 
fluxes of the three components. Sample fits of the [Fe n] profile are 
shown in Paper I. 


3.2 Emission-line flux distributions 

In order to obtain the flux distributions, we fitted only one Gaussian 
to all emission lines, for simplicity. Although two or three compo¬ 
nents reproduce better the profiles at some locations, as discussed 
above, the resulting error in the fluxes at these locations is lower 
than 5%. 

Figure [3] shows the flux distributions for all measured emis¬ 
sion lines. In each panel, black represents masked regions where 
the flux values are smaller than the standard deviation of the con¬ 
tinuum values near the emission line or the uncertainty in the flux is 
higher than 40 %. Light gray contours delimit the masked regions 
and the gree n contours overlaid to som e panels are from the ra¬ 
dio image of lUlvestad & Wilsonl dl989h . At most locations the un¬ 
certainties in flux are smaller than 15 %. For illustration purpose, 
we have rebbined the spaxels to 1/3 of their original size and then 
interpolated their fluxes. As the original spaxels (0705 x0705) are 
smaller than the angular resolution of our observations, this proce¬ 
dure does not affect the spatial resolution of the maps significantly. 
The line emission is more extended along the PA = 60/240°, ex¬ 
tending up to 175 to both sides of the nucleus, while to the perpen¬ 
dicular direction (PA = -30/150°) the emission is extended to 077 
from the nucleus. The flux distributions of all emission lines show 
a good correlation with the radio structures, with the two peaks 
of emission associated to the southwestern and northeastern radio 
hotspots. Some differences are observed among distinct flux dis¬ 
tributions. While the [Fe n] and H 2 fluxes peak at the location of 
the southwestern radio structure at 076 from the nucleus, the Hi 
recombination lines present their highest fluxes at the location of 
the northeastern radio hotspot at 0750 from the nucleus. Another 
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Table 1. Measured emission-line fluxes for the five positions along PA = 60° and PA = -30° within 0735x0735 aperture in units of 10 16 ergs 1 cm 2 . 


dvac (A) 

ID 


075ONE 

0760SW 

Nucleus 

0735SE 

0735NW 

11471.3 

[P II] 1 D 3 - 3 P i 


12.1+5.3 

6.3+4.3 

- 

- 

- 

11665.6 

[Fen] b 4 P\/2~ 

a 2 P3/2 

21.1+4.8 

- 

- 

- 

- 

11863.7 

[Fen]6 4 Z> 1/2 - 

b A F 5 , 2 

10.4+5.7 

- 

- 

- 

- 

11886.1 

[P II] l D 2 - 3 Pi 

18.3+4.4 

26.9+8.6 


- 

- 

12570.2 

[Fen] a 4 Di/2 - 

a b Di :) /2 

179.4+15.3* 

299.9+18.0** 

1 10.2+8.2 

35.0+8.5* 

110.1+9.3* 

12791.2 

[Fen] a 4 D 3/2 - 

a 6 D 3/2 

9.9+2.3 

4.7+1.2 

- 

- 

14.6+9.5* 

12821.6 

HiPa 


117.3+5.2 

106.1 + 10.8 

79.0+9.3* 

29.1+4.4* 

157.9+12.1* 

12946.2 

[Fen] a 4 D 5/2 - 

a 6 D 5/2 

13.4+4.4 

5.5+3.3 

18.6+9.5* 

0.70+0.13* 

17.8+14.1* 

13209.2 

[Fen] a 4 Di/2 - 

a 6 D 2 /2 

60.9+20.4* 

69.0+7.8 

13.6+8.2* 

2.21+0.26* 

24.5+20.9* 

13281.4 

[Fen] a 4 D 5/2 - 

a 6 D 3/2 

18.8+11.2 

- 

- 

- 

- 

21218.3 

H 2 2 1-0S(1) 


18.7+ 1.7 

45.7+ 1.5 

27.1+7.4* 

16.8+4.8* 

16.5+3.0* 

21542.0 

H 2 1-0 S(2) 


1.9+ 0.7 

2.1+0.6 

- 

3.8+2.9* 

- 

21661.2 

HiBry 


16.4+ 1.5 

14.5+ 1.3 

12.9+2.4* 

10.5+3.0* 

12.5+3.5* 

22234.4 

H 2 1-0 S(0) 


9.4+ 0.7 

11.4+ 0.6 

- 

- 

- 

22477.6 

H 2 2-1 S(l) 


4.9+0.62 

5.2+ 1.0 

- 

- 

- 

24084.7 

H 2 1-0Q(1) 


21.8+ 1.3 

39.5+ 0.7 

25.8+6.8* 

13.8+6.4* 

15.5+6.0* 

24136.7 

H 2 1-0 Q(2) 


6.3+2.5 

7.9+ 1.9 

- 

- 

- 

24218.0 

H 2 1-0 Q(3) 


18.0+ 2.8 

32.3± 2.8 

25.0+14.4* 

15.4+9.9* 

17.6+13.1* 

24369.7 

H 2 1-0 Q(4) 


- 

12.6+ 1.9 

- 

- 

- 

24548.5 

H 2 1-0 Q(5) 


14.3± 3.6 

22.9+ 4.7 

- 

- 

- 

24755.5 

H 2 1-0 Q(6) 


- 

8.4+ 4.1 

- 

- 

- 

25000.7 

H 2 1-0 Q(7) 


- 

10.5+8.2 

- 

- 

- 


* Emission-line profile fitted by two Gaussian curves and the flux represents the sum of the fluxes of both components; 
** Emission-line profile fitted by three Gaussian curves and the flux represents the sum of the fluxes of the components. 


difference is that the H 2 emission is less collimated and more ex¬ 
tended in all directions, as most clearly seen in the H 2 >12.1218 flux 
distribution (that presents the highest signal-to-noise ratio among 
the H 2 lines). 


3.3 Line-ratio maps 

The excitation mechanisms of the [Fe n] and H 2 emission 
lines can be investigated using emission-line ratio maps. Fig¬ 
ure 0| shows the [Fen]/ll.2570yi/m/[Pn]/ll.8861 /dm (top-left 
panel), [Fe n]/l 1.2570//m/Pa (3 (top-right panel), Fl 2 /12.1218 yum/Bry 
(bottom-left panel) and Fl 2 T2.2447yum/F[ 2 T2.1218yi/m (bottom- 
right panel) line ratio maps. The first two are useful to investigate 
the [Fen] emission origin and the last ones to investigate the H 2 
emission origin. 

The [Fe n]/[P n] map shows values ranging from ~2 to 10, with 
the highest values observed at the location of the northeastern radio 
knot and the lowest values at regions next to the nucleus of the 
galaxy. At the position of the southwestern radio knot, [Fe n]/[P n] 
is ~7. The [Fen]/Pa/? line ratio shows values ranging from 0.5 to 
up to 3.5, with the highest values observed in regions next to the 
southwestern radio knot. A small enhancement in this ratio is also 
observed next to the northeastern radio knot, where the values are 
approximately 2.5. At most locations, the [Fen]/Pa/3 ratio shows 
values in the range 0.6 to 2.0. Values smaller than 0.6 are observed 
only in a small region at 075 north-northeast of the nucleus. 

The H 2 /Bry line ratio shows values in the range from 0.2 to 
4.5, with the highest values observed in regions surrounding the 
southwestern radio knot and at the nucleus. For most locations, the 
values of H 2 /Bry are in the range 0.6-2.0 and the smallest values 
of 0.3 are observed at 075 north-northeast of the nucleus, at the 
same location where the [Fe n]/Pa/? ratio map presents its smallest 
value. Finally, the H 2 T2.24/2.12 ratio map presents values ranging 


E(B-V) 

fi i i | i i i i | i i r 



1 0 -1 


arcsec 


2.00 

1.50 

1.00 

0.50 

0.00 


Figure 5. Reddening map obtained from the Pa/3/Bry line ratio. Green con¬ 
tours are from the radio image and the central cross marks the position of 
the nucleus. 


from 0.1 to 0.3 with the highest values seen at 075 northeast of the 
nucleus, associated to the radio knot there. 

In Figure [5] we present a reddening map obtained from the 
Pa/?/Bry line ratio using 


E(B-V) = 4.74 log 


5.88 


Fpap/FBry / 


i 


( 1 ) 


where F Pa/3 and F Bry are the fluxes of Paj3 and Bry emission 
lines, respectively. We have adopted the intrin sic ratio F PaP IF Br y = 
5.88 corresponding to case B recombination lOsterbrock & Ferlandl 
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Figure 3. Integrated flux maps for the emission lines identified in the top of each panel, in units of 10“ 17 erg s -1 cm -2 . The central c ross m arks the position of 
the continuum peak. The green contours overlaid on some maps are from the 6 cm radio continuum image from luivestad & Wilsonl 11989). 


[2006) and used the reddening law of ICardelli. Clavton & Mathis! 
dl989h . The values of E(B - V ) for NGC5929 range from 0 to 2, 
with the highest values observed to the south, southeast and east of 
the nucleus, while smaller values are observed approximately along 


the region covered bye the radio jet at PA ~ 45°, from the northeast 
to the southwest of the nucleus. 
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3.4 Stellar Kinematics 

The K-band spectra of Fig. [T] show clearly the CO absorption band 
heads around 2.3 yi/m. We used the Penalised Pixel-Fitting (ppxf) 
method of lCappellari & Emsellem (2004 ) to fit these bands and ob¬ 
tain the line-of-sight velocity distributions (LOSYD) of the stars, 
using the Gemini library of late spectral type stars observed with 
the Gemini Near-Infrared Spe ctrograph (GNIRS) IFU and NIFS 
dWinge. Riffel & Storchi-Bergamann|[2009 ) as templates. The stel¬ 
lar LOSVD was approximated by a Gaussian distribution. The ppxf 
outputs the stellar radial velocity (V*), the corresponding velocity 
dispersion (cr*), as well as the uncertainties for both parameters at 
each spaxel. 

Figure[6]presents the resulting maps for V* and cr*. The signal- 
to-noise ratio of the spectra at locations next to the borders of the 
field of view was too low and we could not obtain good fits. These 
borders were then masked out. The white (black) regions in the V* 
(cr*) map correspond to these locations, where the uncertainties in 
V* or cr* are larger than 30 km s -1 . The stellar velocity field shows a 
velocity amplitude of about 100 km s -1 , with redshifts to the north¬ 
east and blueshifts to the southwest, thus opposite from the rota¬ 


tion field observed for the gas, that has redshifts to the southwest 
and blueshifts to the northeast, as observed in Fig. [7] The cr* map 
shows values ranging from 40 to 180 km s -1 with a median value of 
cr* = 114 + 7 km s -1 and presents a partial ring of lower cr* values 
(60-100 km s -1 ) with radius of 075 surrounding the nucleus. 


3.5 Gas Kinematics 

We used the centroid wavelength of the emission lines 
[Pn] 21.1886yt/m, [Fen] 21.2570yum, Pa/? and H 2 T2.1218yt/m at 
each position to map the velocities of the ionized gas forbid¬ 
den lines, ionized gas permitted lines and molecular gas. These 
lines have been chosen because they present the highest S/N ratios 
among their species. 

These maps are shown in Figure [7j where the white regions 
have been masked due to bad fits, following the same criteria used 
for the flux distributions of Fig. [3] All velocity maps show redshifts 
of up to 220kms _1 to the southwest of the nucleus and blueshifts 
of similar amplitude to the northeast of it. The zero velocity line 
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galaxy spectra. 




105. 


180. 
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is oriented approximately along PA = -30°/150°, coincident with 
the strip where the lines are double. 

Figure [8] presents the corresponding velocity dispersion (cr) 
maps. The cr values for all emission lines are in the range 30 to 
200kms _1 . At most locations, cr is lower than 100 km s -1 . The 
highest cr values are observed along the SE-NW strip (that has 
a width of « 073, thus 50 pc at the galaxy) crossing the nucleus 
perpendicularly to the radio jet and are due to the double compo¬ 
nents observed in the emission lines at these locations. These dou¬ 
ble components were attributed to and equatorial outflow in Paper 
I, where we presented the results only for [Fen]. Now we compare 
the [Fe n] and [P n] maps with those of Pa/? and H 2 and we find 
some differences. The [Fe n] and [P n] maps show that, besides the 
high cr values observed along the SE-NW strip, there is also an in¬ 
crease of cr at and around the locations of the radio hotspots. The 
Pa/? and H 2 maps, on the other hand, do not show any evidence 
of increased cr there. A comparison between the profiles of these 
lines at the location of the radio hotspots can be seen in Fig. [21 it is 
clear that the profile of the [Fe n] line is broader than those of Pa/? 
and H 2 lines as already reported in Paper I, being attributed to the 
interaction of the radio jet with the ambient gas. 


3.6 Channel Maps 

In figures [9] [TO] and QT] we show velocity-channel maps along 
the Pa/?, [Fe n] and H 2 emission-line profiles, respectively. The H 2 
channel maps show emission for a velocity range from ~ -300 to 
300 km s -1 , with the emission moving from northeast to southwest 
as the velocity increases, being consistent with a rotating disk. At 
low velocities, the Pa/? channel maps are similar to that observed for 
H 2 , but additional Pa/? emission from high velocity gas is observed 
associated to the radio hotspots to the northeast at the blueshifted 
channels and to the southwest at at the redshift channels. This sug¬ 
gests an interaction of the radio with the Pa/? emitting gas. For the 
[Fe n] emission, the correlation of the line emission with the radio 
hotspots is more clear with emission being observed at velocities 
of up to -560 km s -1 at blueshifts and 500 km s -1 at redshifts. 


4 DISCUSSION 


4.1 The emission-line flux distributions 


All emission lines present their flux distributions more extended 
along the PA - 60/240°, which is the or ientation of the radio jet 
(lUlvestad & Wilsonll 19891 : ISu et al.lll996h . The companion galaxy 
NGC 5930 is also at PA = 60. Along the perpendicular direction, 
the line-emission extends only to ~077 from the nucleus for most 
lines (a bit more for H 2 ). In particular, along the SE-NW strip, we 
have concluded in Paper I that the gas is outflowing from the nu¬ 
cleus. 

Although the bulk of the emission-line flux distributions is 
similar for all species, some differences are observed: (i) while 
the flux in the Hi recombination lines peak at 075 northeast of 
the nucleus, the [Fen], [Pn] and H 2 lines peak at 076 south¬ 
west of the nucleus; (ii) the [Fen], [Pn] and Hi lines are more 
collimated along the radio axis than the H 2 emission. We have 
similar differences between the flux distribution (and kinemat¬ 
ics) of the ionized and molecular gas for most of the Seyfert 


galaxies we have studied so fa r in this project (Riffel et al. _ 

2008, 20091, 2014al: Riffel. Storchi-Bergmann & Nagad I 2 OIO 1 


Riffel & Storchi-Bergmann _ 2011 alibi: IStorchi-Bergmann et alj 

2009, 2010; Schonel et al. 2014). 


We also note that the flux values along the SE-NW strip are 
lower than along the direction of the radio jet, revealing that the out¬ 
flowing gas discovered in Paper I does not contributes much to the 
total luminosity. This is consistent with the interpretation that the 
radiation from the AGN escapes mainly along the ionization axis, 
which seems to be coincident with the radio axis and is blocked 
along the perpendicular direction due to the dusty torus postulated 
by the AGN unified model. This interpretation is supported by our 
discussion in Sec. 14.71 in which we conclude that X-ray heating (as 
X-rays can escape through the torus) is the main excitation mech¬ 
anism of the H 2 emission of both kinematic components observed 
along PA = -30/150°. 

We can compare our results with previous ones from the lit¬ 
erature. Optical Int egral Field Spectroscopy (IFS) was obtained 
for NGC 5929 by iFerruit et aD dl997h at a seeing of 078 and by 
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Figure 7. Centroid velocity maps derived from th fit of the [Pn] dl.l886yum, [Fen] 41.2570yum, Pa J3 and H 2 42.1218//m emission-line profiles using a single 
Gaussian. The central cross marks the position of the nucleus and the scale of the color bar is shown in units of km s -1 . 


IStoklasova et alJ d2009h with the instrument OASIS at the Canada 
France Hawaii Telescope ( CFHT) in the low spatial resolution 
mode (074 sampling). iFerruit etldl dl997h present a contour map 
for the [N n]+Hor emission, showing that it is more extended along 
the PA of the radio jet and with the peak of emission observed at 
~076 southwest of the nuc leus, in good agreement with our flux dis¬ 
tributions in the near-IR. IStoklasova et al.1 d2009h present flux and 
kinematic maps for the [O iii]/15007, H/3 and fN i]TT5198, 5200, 
finding a similar flux distribution to that of lFerruit et alJ (1 19971 ). Al¬ 
though the optical emission-line flux distributions are also more ex¬ 
tended along PA = 60/240°, similarly to what we observe, the opti¬ 
cal lines do not show an emission peak to the northeast, as we have 
observed for the near-IR lines. We attribute this differences to the 
worse spatial resolution of the opticcdobservations. This interpreta¬ 
tion is supported by the work of lRosario et al.1 (2010), who present 
a Hubble Space Telescope narrow-band image for the [Om ]+HJ3 
emission that clearly shows two main emission regions, one to the 
northeast and other to the southwest, with the highest fluxes ob¬ 
served to the southwest, in good agreement with the [Fen], [Pn] 
and H 2 flux distributions. 

Why do the near-IR Hi recombination lines present dis¬ 
tinct flux distributions from that of H p in the flux map of 


IStoklasova et al. ( 2009 )? As seen in Fig. [3] the Pa fi and Bry flux 
maps present their peak at 075 northeast of the nucleus with a sec¬ 
ondary peak at 076 southwest of it. On the other hand, the H/? flux 
distribution shows its emission peak to the southwest and although 
the emission is also extended to the northeast, there is no emis¬ 
sion peak as observed for the near-IR lines. We attribute this dif¬ 
ference to dust extinction, that is stronger in the optical than in 
the near-IR. Indeed, as seen in Fig. [5j the reddening is somewhat 
larger to the east/northeast of the nucleus supporting the conclu¬ 
sion that the near-IR emission is probing a dustier region than the 
optical emission. A difference of 0.5 mag in E(B - V ) between 
the northeast and southwest sides corresponds to a ratio of 2.2 
b etwe e n the Pa 0 flu xes o f the northeast and southwest using the 
ICardelli. Clavton & Mathis ( 1989 ) extinction law, while for HJ3 the 
ratio would be higher, of « 2.9. This difference would be enough 
to make the northeastern hot spot brighter than the s outhwestern 
hotspot in the H ft flux map of Stok lasova et al.1 (120091) . supporting 
our conclusion that the difference between the near-IR and optical 
flux distributions is due to extinction. 
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Figure 8. Velocity dispersion maps derived from the [Pn] /U.1886jum, [Fen] /ll.2570//m, Pa/3 and H 2 d2.1218//m emission line profiles. The central cross 
marks the position of the nucleus, the scale of the color bar is in units of km s -1 and the green contours are from the 6 cm radio image of luivestad & Wilsonl 
dl989h. 


4.2 The stellar kinematics and the mass of the SMBH 


The stellar velocity field (left panel of Fig. n shows a rotation 
component, with the north-eastern side of the disk approaching 
and the south-western side receding with a velocity amplitude of 
about 80 kms -1 . In order to derive the systemic velocity ( V s ), ori¬ 
entation of the line of nodes ('Fq), inclination (/) and eccentricity 
(e) of th e disk, we fitted the stellar velocity field using the disk- 
fit code (Spekkens & Sellwoodll2007l : ISellwood & Sanche3l2010l : 


fit code (bpekkem 

iKuzio et al.l2012l) . 


= - <0 500 * 0 049 > + <4 42010 295)108 (so^f) ■ 

We estimated a mass of M. = 5.2^-f 2 x 10 7 M©. 

The central velocity dispersion quoted in the Hyperleda 
database ( Paturel et a l. 2003 ) is cr* = 120.6 + 12.9 kms -1 , and re¬ 
sults in M # = 3.4^ 9 x 10 7 M 0 . 


4.3 The gas Kinematics 


The best rotating disk model is shown at the top-central panel 
of Fig. [13] The observed velocity field is shown at the top-left panel 
and the residual map is presented at the top-right panel of the same 
figure, respectively. The corresponding kinematic parameters are 
shown in Table |2] 

The stellar velocity dispersion (cr*) can be used to estimate 
the mass of the supermassive black hole (M.) at the center of 
NGC5929, using the M. - cr relationship. We obtain cr* = 133 ± 
8 kms -1 using the pPXF for an integrated K spectrum over the 
who le field of view. The M. can be estimated by (iKormendv & Hoi 
l2013h : 


4.3.1 Equatorial outflow 

In Paper I we have reported the discovery of a peculiar gas outflow 
along the SE-NW strip that has a width of « 50 pc and an extent 
of ~300pc along PA = -30/150°, perpendicularly to the radio 
jet. At locations away from the SE-NW strip, the gas kinematics 
is consistent with orbital motion in a disk that is counter-rotating 
relative to the stellar disk. This rotation component is clear in our 
Fig. |7] that is based on a single-Gaussian fit to the emission line 
profiles and can be compared to the stellar velocity field shown in 
the left-panel of Fig. [6] Our stellar and g as vel ocity fields are very 
similar to those presented bv lStoklasova et al l (120091) . and derived 
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Figure 9. Channel maps along the Pa fl line profile for a velocity bin of 75 km s 1 , corresponding to three spectral pixels. The corresponding velocity is shown 
in the top-right corner of each panel in units of km s -1 and the green contours are from the radio jet. The color bar shows the flux scale in logarithmic units. 


from optical IFS, although they did not detect the outflow in the 
SW-NW strip due to their poorer spatial resolution. 

In order to better show the distinct kinematic components 
along the SE-NW strip, where the mass outflow was observed in 
Paper I, we extracted one-dimensional cuts from the [Fe n] flux dis¬ 
tribution, velocity field and cr map for the two components along 
the strip. The cr of the two components was kept the same, as done 
in Paper I. These cuts are shown in Fig.QJ] where it can be seen that 
the velocity of the two components are almost constant along the 
strip, with values of -150 and 150 km s -1 relative to the systemic 
velocity of the galaxy for the blue and red component, respectively. 
The fluxes of both components are similar at most locations, with 
the red component being brighter at distances smaller than 075 from 
the nucleus to the southeast. These components were attributed to 
the presence of an equatorial outflow (perpendicular to the radio 
jet) in Paper I. Such gas outflows appear in recent theoretical mod¬ 


els of accretion disk winds (ILL Ostriker & S unyaevj 


in outflowing torus models fe.g lHonig et al J20li 


2013) as well 


ElitzuJ2pT2l) . 


in which the outflows are originated as a consequence of the con¬ 
servation of the gas angular momentum. 

We can use the observed velocity and geometry of the emitting 
gas to estimate its mass outflow rate (M) in the equatorial outflow. 
The double components observed along PA = -30/150° are at¬ 
tributed to the equatorial outflow from the nucleus. As the outflow 
is observed as two similar velocity components, one in redshift and 


the other in blueshift, we have considered a scenario in which the 
outflow has generated a hollow cylinder of outflowing mass. Con¬ 
sidering the width of the SE-NW strip (50 pc) as being the height 
(h) of the cylinder considered to have a radius R 0 , the mass outflow 
rate can be estimated by 

M = m p N e vfA , (2) 


where m p is the proton mass, N e the electron density, v is the ve¬ 
locity of the outflowing gas, / is the filling factor and A = 2 nR^h 
is the lateral area of the cylinder. As the velocity of the outflow is 
approximately constant along the strip and there is no decreasing 
velocity as a function of distance from the center to NW and SW, 
we conclude that this cylinder has a large radius, at least the size of 
the region where emission is observed, R 0 >075. 

Assuming N e = 500 cm -3 , / = 0.01 (which are typical val¬ 
ues es ti mated for other Seyfert galaxies - e.g. ISchnorr-Muller al.1 
l2014al lbL l201ll : [Storchi-Bergmann et al] l2010l) . v = llOkms -1 
(from Fig. [12]) and i? 0 =075=87.5 p C , we obtain a lower limit 
(as the radius for the cylinder can be considered a lower 
limit) to the mass outflow rate in ionized gas of M > 
0.38 M©yr -1 . This value is in the range of mass outflow rates 
observed for other Seyfert galaxies (Crenshaw & Kraemer ]2007 ; 


Barbosa et al] 120091: IStorchi-Bergmann et al. ~ l2010l : iRiffel et alJ 
2009; Riffel & Storchi-Bergmann 2011a) with similar luminosities, 


although it is much smaller than those found in recent studies for 
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Figure 10. Same as Fig.[9]for the [Feii]/U.2570 /urn emission line. 


high luminosity Seyfert 2 galaxies (log(L[OIII]5007/L Q ) > 8.7) 
which show M out « 370 - 2700 M© yr -1 (e.g. lMcElrov et al ji2014h . 

Finally, due to the apparent unique velocity of the outflow, it 
looks like that it is not continuous, but was generated by a ’’blast” 
that has produced an expanding (cylindrical) shell of gas. If this 
is the case, we can estimate the age of the “blast” by r = R 0 /v 
resulting in r = 0.8 Myr. Actually, this can be considered an upper 
limit for the age, as the outflow may have decelerated since the 
ejection from the AGN. 

4.3.2 Interaction between the radio jet and emitting gas 

Our data allowed us to observe also the kinematic effects of the 
interaction of the radio jet with the [Fen] emitting gas, as discussed 
in Paper I based on a two Gaussian component fit to the [Fe n] line 


profile. The interaction of the radio jet with the gas produces wings 
in the [Fe n] line profile at locations around the radio hotspots. In 
the single-Gaussian fit, this effect appears as an enhancement of cr 
at these locations, as seen in Fig. [8] iRosario etal . ( 12010 ) claims 
also to have found similar signatures of shocks due to the radio jet, 
using optical long slit observations with the HST. 

Although signatures of the interaction of the radio jet with 
the ISM are clearly observed for the [Fe n] emission, the H 2 and H 
recombination lines show weaker signatures, only observed in the 
channel maps. This indicates that the [Fen] emission traces a more 
disturbed kinematics than the molecular and atomic hydrogen. 
This can be interpreted as due to shocks destroying dust grains 
in its passage through the gas, releasing the Fe atoms, that are 
then ionized and emit.This conclusion is supported by the reduced 
reddening observed along the ionization and radio jet axis of the 
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Figure 11. Same as Fig.[9]for the H2/12.1218 //m emission line for a velocity bin of 60 km/s, corresponding to two spectral pixels. 


AGN, and also in agreement with previous results we have found 
for other Seyfert galaxies: the molecular gas and the ionized gas 
emitting the [Fe n] lines trace distinct kinematics and flux distribu¬ 
tions with the former more restricted to the plane of the galaxies 
an d the latt er usually assoc ia ted to outflows^ from the nucleus 


2010al; Riff el & Storchi-Bergmann 

2011a 

Riffel, Storchi-Bergmann & Winge 2013]” Riffel et al. 

2014a: 

Storchi-Bergmann et al. 2009, 2010; Barbosa et al. 2014]). 


We can also speculate on the age of the radio jet. Assuming 
the jet velocity of 0.1 c, if the jet were in the plane of the sky with 
a size of O^, the age would be 2850 yr. If the inclination of the 
jet relative to the plane of the sky is 45°, the age would be about 
4000 years and we thus conclude that the radio jet is younger than 
the ’’blast” generating the equatorial outflow. Alternatively, the jet 
is much more inclined toward us and could have been emitted to¬ 
gether with the ’’blast”, taking also in consideration that the velocity 
of the outflow may have been larger when it was generated. 


4.3.3 Rotating disk 

As observed in the stellar velocity field, the gas velocity fields 
(Fig. 0 also support the presence of a rotating disk compo¬ 
nent, although the gas is clearly counter-rotating_relative to the 
stars. This was already observed by IStoklasova et al.1 J2009h . as 


mentioned above. In order to confirm the presence of this com¬ 
ponent and derive the corresponding kinematic parameters, we 
fitted th e Pa J3, fFenl and H 2 velocity fields using the disk- 
fit code (Spekkens & Sellwoo3l2007l : ISellwood & Sanchezll2010l : 
iKuzio et al.ll2012h for a rotating disk model. In order to fit the rota¬ 
tion component of the gas velocity field, we excluded the velocities 
from the SE-NW strip. The exclusion followed the same criteria 
adopted in Paper I, in which regions with cr > 100 kms -1 for the 
[Fen] line along PA = -30/150° were attributed to the equatorial 
outflow. 

The corresponding best models are shown in the central pan¬ 
els of Fig. [13] while the observed velocity fields and the residual 
maps are presented in the left and right panels, respectively. As the 
kinematics for distinct emission lines are similar, we show only the 
maps for the H 2 in this figure. It can be seen that the residuals are 
smaller than 50 km/s at most locations, indicating that the velocity 
fields are reasonably well representing by rotation in a disk. In Ta¬ 
ble [2] we present the kinematic parameters for the best fit models. 
The systemic velocity (V 5 ), ellipticity of the orbits ( e ), inclination 
of the disk (/) and position of the kinematical center (X cen , Y cen ) 
for the molecular and ionized gas and stars are very similar. The 
kinematical center is measured relative to the peak of the K-band 
continuum and was kept fixed for the stellar velocity field in order 
to reduce the number of free parameters to be fitted, as the stellar 
velocity field is nosier than those for the emission lines. 

The orientation of the line of nodes ( v Fo) obtained from the fit 
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SE Position (arcsec) NW 

Figure 12. One-dimensional cut for the fluxes (top), velocity (middle) and 
velocity dispersion (bottom) from two-Gaussian components fit to the [Fe n] 
emission-emission line profile along the SE-NW strip. The velocity disper¬ 
sion was contained to be the same for both components during the fits. 


of the gas velocity fields are consistent with each other, while T'o 
for the stellar velocity field is approximately opposite to that of the 
gas (differing by 190°). 

We can compare the kinematical parameters from Table[2]with 
those for the large scale disk, that are quoted in the Hyperleda 
database (iPaturel et al .11 20031) . The heliocentric systemic velocity 
obtained from the fit is in good agreement with the large scale value 
(V s ~ 2503 kms -1 ) obtained from optical emission lines, while 
the inclination of the disk and the ellipticity of the orbits are a bit 
smaller than those for the large scale disk (/ « 24°, e - (l -b/a) ^ 
0.09). The orientation of the line of nodes is 10 - 20° smaller than 
that of the large scale disk (e.g. lSchmitt et al.lll997 ). We attribute 
the origin of this rotating gas disk to the interaction with NGC5930. 
The gas velocity amplitude, corrected by the inclination of the disk 
is about 580 km s -1 .The high rotation velocities indicate that this 
gas is still not in orbital equilibrium in the galaxy gravitational po¬ 
tential. This gas is probably the source of the feeding of the AGN. 


4.4 Mass of ionized and molecular gas 

The mass of the ionized gas in the inner 520x520 pc 2 can be esti¬ 
mated by 


M hu w 3 x 10 17 


FBry 

erg s -1 cm -2 



[Mo], 


( 3 ) 


where F Brr is the integrated flux for the Bry emission line 
and d i s the distance to_NGC5929 (lOsterbrock & Fcrland 2006 : 
IStorchi-Bergmann et al.l l2009). We have assumed an electron tem¬ 
perature T = 10 4 K and electron density N e - 10 2 cm -3 . 

The mass of the warm molecular gas is given by 


where F Hl x 2.1218 is the integrated flux for the H 2 T2.1218jum emis¬ 
sion l ine and we have used the vibrational temperature T=2000K 
(e.g. lScoville et alll982l : lRiffel et al.ll2008h . 

Integrating over the whole IFU field (520x520 pc 2 ) we obtain 
F Brr = 3.4 + 0.4 x 10 -15 erg s -1 cm -2 and F H2 a 2.1218 ~ 7.3 + 0.9 x 
10 -15 ergs -1 cm -2 , that results in M HII - 1.3 + 0.2 x 10 6 M© and 
M H2 = 471 + 58 M©. 


The masses of warm molecular and ionized gas are in 
the range of values found for other active galaxies, which 


have Mi 


HII 


= 10 4 


10 7 M© and M, 


ih- 


= 10 1 - 10 3 M© 


(Riffel et ah _ 20081. 120091: Riffel. Storchi-Bergmann & N agar 

2010al: IRi ffel & Storchi-Bergmann 1201 lbl : IRiffel et al. I l2014aT 


Storchi-Bergmann et al .1 [2 009: Riff el. Storchi- Bergman n & Wingei 
20131 : iMazzalav et al.1120131) However, the galaxies from these 


studies are located at a large range of distances (10-100 Mpc) 
and thus the corresponding sizes of the regions covered by the 
observations - obtained with NIFS (and SINFONI) - are distinct 
for the distinct objects. 

In order to be able to compare similar quantities, we now 
calculate the average mass surface density, instead of the mass. 
For NGC5929 we get Y,m Hi = 1-7 x 10 -3 M©/pc 2 and Y,m hu = 

4.8 M©/pc 2 for the mass surface density of the warm molecular 
and ionized gas, respectively. 

These values are in the range observed for other galax¬ 
ies, that renge from Tjm h , = 2.4 x 10 -4 M©/pc 2 for NGC 1068 
( [Riffel et al.[l2014al) to = 7.1 X 10 -3 M©/pc 2 for NGC 2110 
(IDiniz et al.ll2015l) for the warm m olecular gas and fr om Yjm hii - 

1.8 x 10 _1 M o /pc 2 for NGC 106 8 feiffel et al.ll2014ah to = 

42.8 M 0 /pc 2 for NGC 4151 dStorchi-Bergmann et al.lbooT for the 
ionized gas. However, as the galaxies of our sample present dis¬ 
tances in the range of 9.3 Mpc (for NGC 4151) to 93.8 Mpc (for 
Mrk79), distinct average densities are expected for similar radial 
density profiles with decreasing gas densities as a function of dis¬ 
tance from the nucleus. 

In order to compare the masses of ionized and and warm 
molecular gas closest to the nucleus, we used the previous data 
obtained by our group to estimate the masses within a fixed phys¬ 
ical aperture of 100 x 100pc 2 . This aperture was chosen as be¬ 
ing the whole field of view of the nearest object of the sample. 
The mass of ionized gas ionized gas within this aperture is in the 
range 0.3 - 27 x 10 5 M© and the mass of warm molecular ranges 
from 3 to 455 M©, approximately. The lowest values for both ion¬ 
ized and molecular gas masses are observed for NGC 4051, while 
NGC 1068 presents the highest values. For NGC 5929 we find 
M H u - 3 + 0.5 X 10 5 M© and M Hl =47 + 6 M© for the inner 
100 x 100pc 2 , which corresponds to an aperture of only 076x076. 
The median ratio between ionized and molecular gas masses is 
6.9 x 10 3 , similar to the ratio for NGC 5929 of ~ 6.4 x 10 3 . 

However, the total mass of molecular gas, including the cold 
gas is much larger than the values obtained here. The ratio be¬ 
tween cold and warm molecular gas observed in the central re¬ 
gion of active galaxies is in the range 10 5 -10 7 (iDale et ahl120051 : 
iMuller-Sanchez et aPl2006l : IMazzalav et al.ll2013l) and thus, the to¬ 

tal amount of molecular gas in the inner 3" x 3" of NGC 5929 
should be at least 4.7 x 10 7 M© with a surface mass density of 
170 M©/pc 2 . 


4.5 The origin of the [Fe n] emission 


M Hi 


* 5.0776 x 10 13 


Fh 2 A 2.1218 
erg s _1 cm -2 



The origin of the [Fen] emission in AGNs can be investigated 
using the line-ratio maps shown in Fig. [4] The [Fe n] emission 
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Figure 13. From left to right: observed velocity field, rotating disk model and residual map (observed - model) for the stellar (top) and H 2 (bottom) velocity 
field. The color bars show the range of velocities in units of kms -1 relative to the systemic velocity of the galaxy and the central cross marks the position of 
the nucleus. 


Table 2. Kinematic parameters from the best rotating disk model. Col 1: Fitted velocity field, col 2: heliocentric systemic velocity, col 3: orientation of the 
line of nodes, col 4: ellipticity, col 5: inclination of the disk and cols 6 and 7: position of the kinematical center relative to the peak of the K-band continuum 
emission. 


Velocity Field 

V s (kms -1 ) 

^0 

e 

i 

(Xcen ? Y ce n) 

h 2 

2500.1+17.4 

25.7°+1.3° 

0.05 + 0.02 

18.4°±1.2° 

(072+071, -071+071) 

Pa/3 

2472.8+18.5 

36.5°±1.1° 

0.05 + 0.02 

18.2°±1.3° 

(-071 +071, 07 0+071) 

[Fen] 

2469.3+18.5 

25.4° ±4.4° 

0.05 + 0.04 

18.2°+7.8° 

(071 +071, 07 0+071) 

Stars 

2489.4+17.1 

219.9°+1.3° 

0.05 + 0.03 

18.2° +6.7° 

(0, 0) fixed 


is excited in partially ionized gas regions, which can be pro¬ 
duced b y X-ray (e.g. Simpson et al.1 Il996h and/or shock heat¬ 
ing (e.g. iForbes & Wardl 1 1993b of the gas. Values of [Fen]/Pa/3> 
2.0 indicate that most of the [Fen] line emission is pro¬ 
duced by shocks, while [Fen]/Pa/3< 0.6 indicate that photo- 
ionization dominates the emission (Rod riguez-Ardila et al.l 120041 : 

iRodrfguez-Ardila. Riff el & PastorizJTo05l) 


The [Fen] T1.2570yum/Pa/? line ratio for NGC5929 is shown 
in the top-right panel of Fig.|4] At most positions this ratio has val¬ 
ues larger than 0.6, with the exception of only a small region at 
~078 north-northeast of the nucleus. At locations co-spatial with 
the radio knots, the [Fen] T1.2570//m/Pa/3ratio is larger than 2.0, 
suggesting that the interaction of the radio jet with the gas has an 
important role in the observed [Fen] emission. At the northeast¬ 
ern radio hotspot, [Fen]/Pa/te 2.3, while at the southwestern ra¬ 
dio hotspot [Fen]/Pa/3^ 2.8. At the nucleus and along THE SE- 


NW strip a typical value for [Fen]/Pa/3 is 1.5, suggesting that both, 
photo-ionization and shocks contribute to the excitation of [Fen], 
with a larger contribution of shocks though (as this value is closer 
to 2.0 than to 0.6). 

The [Fen]/ll.2570//m and [Pn]/I1.8861 /dm lines have simi¬ 
lar excitation temperatures, and their parent ions have similar ion¬ 
ization potentials and radiative recombination coefficients and the 
[Fen]/[Pn] line ratio map shown in top-left panel of Fig.|4]is also 
useful to investigate the [Fe 11 ] emission origin. Values larger than 
2.0 indicate that shocks have passed through the gas destroying 
the dust grains, releasing the Fe and thus enhancing its observed 
abundance. For supernova remnants, where shocks are the dom¬ 
inant e xcitati on mechanism, [Fe n]/[P 11 ] is typically higher than 
20 i lOliva et alJlioOlh . In NGC 5929, at all locations where both 
lines are detected, [Fen]/[Pn] is larger than 2.0, suggesting that 
shocks contribute to the [Fen] emission and supporting the con- 
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elusion reached from the [Fe n]/Pa/3 ratio map. At the position of 
the hotspots, values of up to 10 are observed, supporting an even 
stronger interaction of the radio jet and with the gas. 

As the [Fe n] emission-line profiles are complex at many loca¬ 
tions, as seen in Fig. [2] with extended wings at locations close to 
the radio knots and double components along the SE-NW strip, we 
also constructed [Feii]/Pa/? line-ratio channel maps (constructing 
line ratios at different velocity bins) in order to determine the gas 
excitation mechanism at distinct velocities. The resulting channel 
maps are shown in Fig. [14] Each panel presents the [Fe n]/Pa/? line- 
ratio map centered at the velocity shown in the top-right corner of 
the panel, relative to the systemic velocity of the galaxy. In order 
to avoid spurious features, we masked regions in which one or both 
lines present fluxes smaller than 3 times the standard deviation of 
the continuum next to the line. These regions are shown in gray in 
the figure and the contours are from the radio image. 

The [Fe ii]/Pa/3 ratio shows values ranging from 0.2 to up to 
4.5. The highest values are observed for the gas at the highest ve¬ 
locities (largest blueshifts and redshifts) and present an excellent 
correlation with the radio knots. These high [Fe n]/Pa/? values are 
interpreted as a strong evidence that the [Fe n] is produced by the 
interaction of the radio jet with the NLR gas. For smaller velocities, 
at most locations the values are smaller than 2.0, with typical values 
of ~ 1.2, suggesting that X-rays from the central AGN are the main 
drivers of the [Fe n] emission. Line emission is observed along the 
SE-NW strip only at velocities lower than 180kms _1 . The values 
of [Fe ii]/Pa/5 there are somewhat higher the one observed at other 
locations, but still much smaller than those seen at the radio knots. 
Typical values are [Fen]/Pa/?= 2-3, suggesting that shocks are 
also important along the SE-NW strip. These shocks may be asso¬ 
ciated to the equatorial outflows we have reported at this locations 
in Paper I. 

We conclude that the main excitation mechanism of the 
[Fen] emission in NGC5929 are shocks due to the radio jet, 
with some contribution from X-ray heating at locations distant 
from the radio structures. The highest velocity gas is associated 
to the interaction of the radio jet with the ISM and thus we 
conclude that its [Fen] emission is produced mainly by shocks, 
while X-ray excitation is more important for the rotating gas 
at lower velocities. These results are supported by the line-ratio 
maps and by the correlation between the emission-line flux dis¬ 
tributions and the radio image, and are in good agreement with 
those of previous similar studi es for other Seyfert galaxies (e.g. 

Rodriguez-Ardila et al. 2004; Rodriguez-Ardila, Riffel & Pastoriza 
2005: iRiftcl et al . 2006; Riffel Stor chi-Bergmann & Nagadl2010al : 

Storchi-Bergmann et all 19991 2009h. 


4.6 H 2 excitation 


The excitation of the H 2 line emission in the near-IR can 
be due to (i) fluorescence through absorption of soft-UV 
photons (912-1108 A) in the Lyman and Werner bands 
felack & van Dishoeckl 1987) and/or by thermal processes 
due to the heating of the gas by shocks (due to the interac¬ 
tion of the radio jet with t he emitting gas and/or supernovae) 


( Hollenbach & M cKe e 1989) or by X-rays from the central AGN 


dMalonev. Hollenbach & Tielenslll996l) . Several works have been 


(e-g- 

Reunanen, Kotilainen & Prieto 2002^ Rodriguez-Ardila et al. 

2004 

; Rodriguez-Ardila^,Riffel & Pastoriza 2005; Davies et al. 

2005 

; Ramos Almeida, Perez Garcia & Acosta-Pulido 

2009 

) and integral-field spectroscopy (e.g. Riffel et al. 
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Figure 15. Relation between /V upp = and £ upp = 7j for the H 2 emission 
lines for thermal excitation at the locations of the radio hotspots. Ortho 
(S) transitions are shown as filled circles and para (Q) transitions as open 
circles. 
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In order to distinguish between fluorescence and thermal pro¬ 
cesses, we used the observed fluxes for all H 2 emission lines in the 
K band to calculate the thermal excitation temperature. If the H 2 
emission is dominated by thermal processes and under the assump¬ 
tion of an ortho:para abundance ratio of 3:1, the following relation 
is valid JWilman. Edge & Juhnstonell2005h : 

log (^ ) = constant - -p-, (5) 

\ ^igi) T txc 

where F t is the flux of the i th H 2 line, 4/ is its wavelength, A t is the 
spontaneous emission coefficient, g t is the statistical weight of the 
upper level of the transition, T t is the energy of the level expressed 
as a temperature and T exc is the excitation temperature. Thus, if 
the observed fluxes can be reproduced by the equation above, the 
H 2 may be dominated by thermal processes. The resulting plot for 
A upp = (plus an arbitrary constant) vs £ upp = T t is shown in 
Fig.[l5]for the fluxes shown in Table [l]for the locations of the two 
radio hotspots at 075 northeast and 076 southwest of the nucleus. 
We do not show plots for other positions of Table Q] because only 
few H 2 lines were detected at these locations. As seen in Fig. US] 
the observed fluxes are well reproduced by the equation (shown as a 
continuum line) at both positions, indicating that thermal processes 
are the main excitation mechanism of the H 2 lines. The resulting ex¬ 
citation temperature are T exc = 2583 + 40Kand T exc = 2211 + 37K 
at 075 northeast and 076 southwest of the nucleus, respectively. 

The H 2 42.1218/Bry emission-line ratio map shown in Fig. 0] 
can also be used to investigate the nature of the H 2 exci¬ 
tation. For Starburst galaxies and H n regions this ratio is 
H 2 /Bry<0.6. Larger values are observed for Seyfert nuclei, where 
the heating of the gas by shocks and X-rays provide additional 
thermal excitation. Typical values for Se yfert g alaxies are in 
the range 0.6<H 2 42.1218/Bry<2.0 iRodriguez-Ardila et al . 2004 ; 
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Figure 14. Velocity channel maps for the [Fen ]A 1.2570 fdm/Pafi line ratio. Each panel shows the fluxes centered at the velocity shown at its top-right corner 
relative to the systemic velocity of the galaxy. Grey regions are locations where one or both lines present fluxes smaller that 3 cr (the standard deviation of the 
continuum next to the line) and the green contours are from the radio image of lUlvestad & Wilson! dl989l) . 


[Rodriguez-Ardila. Riff el & Pastorizall2005h . although recent stud¬ 
ies suggest a broader range of values_0.4<F[ 2 /l2.1218/Bry<6.0 for 
AGNs, including LINERs feiffel et alJl2013 T NGC 5929 presents 
values larger than 2.0 at most locations, confirming that thermal 
processes are dominating the H 2 excitation. The only exception 
is a small region with values of ~0.3 at ~078 north-northeast of 
the nucleus, where the [Fe ii]/Pa/? also presents small values. Typi¬ 
cal values to the northeast are H 2 /Bry~1.0, while to the southwest 
values of up to 4.5 are observed in regions surrounding the radio 
knot, suggesting that shocks due to the radio jet contribute to the 
H 2 emission at this location. Similar high values are also observed 
at the nucleus and along the SE-NW strip, where the interaction of 
the outflows seen in Paper I with the ambient gas might enhance 
the H 2 emission. 

We also constructed channel maps for the H 2 /Bry ratio in or¬ 
der to investigate the origin of the H 2 emission at distinct kine¬ 
matics. These channel maps are shown in Fig. M The range of 
velocities is much smaller than those seen for the [Fe n]/Pa/3 chan¬ 
nel maps in Fig. [14] and the maps are much noisier due to a lower 
S/N ratio at the K-band lines. At locations of the radio knots, the 
H 2 /Bry values are larger than 2.0, supporting a the contribution of 
shocks due to the radio jet to the H 2 excitation at these positions, in 
particular at the southwestern radio knot, where the highest values 
of H 2 /Bry are observed at the highest velocities channel maps. 

Finally, the H 2 T2.2477//12.1218 can also be used to dis¬ 
tinguish between excitation of the H 2 by thermal processes 
and fluoresce. For thermal processes, this ratio is ~ 0.1 - 0.2 
and for fluorescent excitation it is ~ 0.55 (e.g. M our il 119941: 


Rcunanen, Kotilaincn & Prieto 20021: iRodrfguez-Ardila et al.1 
2004 IStorchi-Bergmann etal] l2009h . The H 2 d2.2477/d2.1218 


ratio map is shown in Fig. [4] and presents values ranging from 0.1 
to 0.25, supporting the conclusion above that fluorescent excitation 
is not important for NGC 5929. 

Thus, we conclude that the H 2 emission observed in 
NGC 5929 is excited by thermal processes, due to heating of the 
gas by shocks and X-rays from the central AGN. At locations co- 
spatial with the radio jet and along the SE-NW strip shocks play 
an important role as indicated by the enhancement of the line ra¬ 
tios. Away from these locations, X-ray heating may dominate the 
H 2 excitation. 


4.7 The origin of the double components 

Do the double components observed along the SE-NW strip have 
distinct origin than the emission from the disk? We can better in¬ 
vestigate the origin of the emission of the outflowing gas by the 
emission-line ratios. Figure [171 presents the [Fen]/Pa/? and H 2 /Bry 
emission-line ratios for the blue (asterisks) and red (diamonds) 
components along a pseudo slit with 0735 width oriented along the 
SE-NW strip at PA = -30/150°. For the [Fen]/Pa/? both compo¬ 
nents show typical values of Seyfert galaxies, with average values 
of 1.0 and 1.3 for the blue and red components, respectively. The 
H 2 /Bry show some values larger than 2 at some positions, but at 
most locations these ratio show values typical of Seyfert galaxies. 
The average values are 1.3 and 2.1 for the blue and red components, 
respectively. As the line ratio of both components show typical val- 
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Figure 16. Same as Fig.[l4]for the Fl 2 /I2.1218/Bry line ratio. 


PA= — 30° 



-0.5 0.0 0.5 

SE Position (arcsec) NW 


Figure 17. [Fen]/Pa/3 and Fl 2 /Bry emission-line ratios along PA - 
-30/150°. The red component is shown as open diamonds and the blue 
component as asterisks. The dashed line shows the average value for the 
red component and the dotted line represents the average ratio for the blue 
component. 

ues of Seyfert galaxies, the emission of the outflowing gas might 
be due to heating by X-rays from the central AGN. 


5 CONCLUSIONS 

We used Gemini integral field J- and K r band spectroscopy of the 
inner 250 pc radius of the Seyfert 2 galaxy NGC 5929, at spatial 
resolution 20 pc and velocity resolution 40 km s _1 , to map the stel¬ 
lar and gas kinematics and the flux distributions of emission lines 
from the ionized and warm molecular gas. Our main conclusions 
are: 

• The flux distributions for the H 2 , H recombination and forbid¬ 
den emission lines are extended to up to 250 pc, with the highest 
flux levels oriented along PA = 60/240° (the direction of the radio 
jet), and well correlated with the radio emission. The H 2 emission 
is more distributed over the whole field of view, while the emis¬ 


sion of the ionized gas is more collimated along the direction of the 
radio jet. 

• The excitation of the H 2 and [Fen] at all locations is domi¬ 
nated by thermal processes, mainly due to heating of the gas by X- 
rays emitted emitted by the central AGN. Some contribution from 
shocks is observed in the high velocity gas at locations next to the 
radio hotpots, as indicated by the [Fe n]/Pa/? and H 2 /Bry line-ratio 
channel maps. 

• The stellar velocity field shows rotation with an amplitude 
of up to ~200kms -1 when corrected for the inferred inclination 
of 18.3°, and a velocity dispersion reaching 180kms _1 . From the 
M. -cr relationship, we estimated a mass for the supermassive black 
hole of M. = 5.2if 2 x 10 7 M 0 using the median stellar velocity dis¬ 
persion of <x* = 133 + 8 kms -1 . 

• The gas kinematics present three components: (i) a counter¬ 
rotating disk (relative to the stellar velocity field); (ii) an equatorial 
outflow (perpendicular to the radio jet), thus in the plane of the 
torus, with a mass-outflow rate of M > 0.4 M 0 yr _1 and (iii) kine¬ 
matic disturbances observed in association with the radio hot spots 
in blueshift to the northeast and redshift to the southwest, what sup¬ 
ports that the radio jet is tilted towards us at the northeast and away 
from us at the southwest. 

• From the Bry and H 2 emission-line fluxes, we calculate the 
mass of ionized and warm molecular gas of M HU - 1.3 + 0.2 x 
10 6 M© and M Hl = 471 +58 M 0 , respectively. These values corre¬ 
spond to mass surface densities of Tjm H2 = 1.7 x 10 _3 M©/pc 2 and 
Yam hii ~ 4.8M©/pc 2 , which are in the range of values observed for 
other galaxies. 

Our favored scenario for this galaxy is that the interaction with 
NGC 5930 has sent gas towards the nucleus of NGC5929, trigger¬ 
ing the nuclear activity. If the observed equatorial outflow can be 
considered the first blast of the AGN, we estimate that this hap¬ 
pened less than 1 Myr ago. 
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